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ABSTRACT 


When  a  metallic  antenna  is  "printed"  on  the  surface  of  a  dielectric 
sheet,  the  original  radiation  pattern  is  altered  by  the  presence  of  the 
dielectric.   This  report  considers  the  case  of  a  half-wave  dipole  antenna 
mounted  on  the  surface  of  a  thin,  rectangular  dielectric  sheet.   For  this 
case,  the  pattern  changes  due  to  the  dielectric  are  found  to  be  quite 
large,  even  for  very  thin  sheets.   Theoretical  radiation  pattern  calcula- 
tions are  presented  and  they  agree  with  the  experimental  patterns  for  thin 
sheets.   The  results  obtained  provide  a  basis  for  estimating  the  radia- 
tion pattern  perturbations  due  to  the  dielectric  for  a  fairly  wide  range 
of  dielectric  constant  and  sheet  dimensions. 
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1 .   INTRODUCTION 

The  University  of  Illinois  Antenna  Laboratory  some  time  ago  set  up 
a  facility  for  etching  '"printed  circuit"  antennas  on  copper-clad  dielec- 
tric sheets.   Accordingly,  it  was  decided  to  investigate  the  effects  of 
a  finite  dielectric  sheet  when  used  as  the  base  for  a  printed  antenna. 
This  report  is  concerned  with  one  aspect  cf  the  problem--that  of  determin- 
ing the  radiation  pattern  of  a  half-wave  dipole  antenna  which  is  on  the 
surface  of  a  thin,  rectangular  dielectric  sheet  and  located  midway  between 
the  edges  as  shown  in  Figure  1 . 

In  order  to  make  the  experimental  model  more  suitable  for  theoreti- 
cal study,  it  was  decided  to  perform  the  experiments  with  the  dipole 
embedded  inside  the  dielectric  sheet,  midway  between  its  two  faces  as 
in  Figure  2.   This  approach  was  chosen  after  experiments  showed  that, 
for  thin  sheets,  there  is  negligible  difference  between  patterns  taken 
with  the  dipole  inside  the  dielectric  and  those  taken  with  the  dipole 
on  its  surface.   A  further  simplification  resulted  from  an  application 
of  the  theory  of  images  which  permits  a  conducting  plane  to  be  passed 
through  the  center  of  the  dipole  and  perpendicular  to  it,   This  made  it 
possible  for  the  radiation  pattern  measurements  to  be  carried  out  on  a 
ground  plane  as  shown  in  Figure  3., 

So  far,  we  have  decided  on  an  antenna  model  suitable  for  experimental 
purposes.   However,  some  further  simplifications  must  be  made  if  we  are 
to  have  an  antenna  model  suitable  for  theoretical  calculations.   Since 
experiments  showed  that  the  pattern  perturbations  in  the  (x,  y)  plane 
were  much  greater  than  those  in  the  other  two  principal  planes,  it  was 
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decided  to  consider  only  the  (x,  y)  plane  pattern  in  this  report.   This 
not  only  simplifies  the  experimental  work  but  also  suggests  reduction 
of  the  three-dimensional  theoretical  problem  to  a  two-dimensional  one. 
This  can  be  done  by  considering  an  equivalent  mathematical  model  having 
an  infinite  line  source  embedded  in  an  infinite  dielectric  strip  as  shown 
in  Figure  4.   Assuming  the  applicability  of  this  mathematical  model  im- 
plies that  the  dimension  h  (of  the  experimental  model)  must  be  great 
enough  so  that  changing  it  appreciably  has  no  effect  on  the  (x,  y)  plane 
pattern  (see  discussion  in  Chapter  2).   Using  the  theoretical  model 
described,  it  is  possible  to  derive  an  integral  equation  for  the  electric 
field  at  any  point  in  space.   In  this  equation,  the  integration  is  over 
the  cross-sectional  area  of  the  dielectric  strip.   However,  the  number  of 
variables  of  integration  can  be  reduced  from  two  to  one  by  assuming  that 
the  dielectric  sheet  is  very  thin  and  that  there  is  no  field  variation 
across  its  narrow  dimension. 

The  radiation  field  was  calculated  by  using  the  assumption  that  the 
field  in  the  dielectric  sheet  is  the  same  as  the  incident  field  in  free 
space.   Such  an  approximation  permitted  the  integration  to  be  carried 
out  numerically  on  a  digital  computer.   This  approach  is  equivalent  to  using 
the  first  term  in  the  iteration  series  solution  for  the  field  inside  the 
sheet.   Therefore  this  solution  will  be  referred  to  as  the  "iteration 
series'"  solution.   It  is  discussed  further  in  Chapter  3  and  Appendix  A. 

A  consideration  of  the  radiation  patterns  presented  in  Chapter  2 
shows  that  even  very  thin  dielectric  sheets  (about  a  hundredth  of  a  wave- 
length thick)  can  have  a  marked  effect  on  the  dipole  pattern.   Fortunately 
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this  effect  is  predictable  when  the  dielectric  is^thin  as  the  theoretical 

results  show,   At  this  point  some  comment  on  the  word  "thin"  seems  to  be 
appropriate;  a  "thin"  polystyrene  sheet  will  be  understood  to  mean  a  sheet 
less  than  three  hundredths  of  a  wavelength  thick.   As  might  be  expected, 
the  usefulness  of  the  theory  is  limited  not  only  by  the  sheet  thickness 
but  also  by  its  dielectric  constant.   Relative  dielectric  constants  above 
four  or  five  severely  restrict  the  applicability  of  the  theoretical  approach 
described  in  this  report. 

The  experimental  results  shown  are  for  polystyrene  with  a  dielectric 
constant  of  2.55  and  most  of  the  measurements  were  carried  out  at  5000  mc . 

However,  the  graphs  shown  can  be  extended  to  other  parameter  values  by 

i 
use  of  the  fact  that  for  constant  sheet  length  -r  and  constant  angle  <p} 

the  pattern  perturbation  magnitude  is  a  function  of  the  factor  (€  -l)  -r 

where  -r  is  the  sheet  thickness  in  wavelengths  and  €   is  the  relative  di- 
rt, r 

electric  constant.   This  means  that,  for  a  fixed  frequency,  the  same  radi- 
ation pattern  will  apply  to  the  two  cases  €   =  2.5,  t  =  .02X  and  €,      =  4.0, 
t  =  .  0l\.   The  theory  breaks  down  for  thick  sheets  but  the  experimental 
results  indicate  the  pattern  shapes  to  be  expected. 

Another  approach  to  the  solution  of  the  integral  equation  is  helpful 
in  interpreting  the  phenomena  observed.   The  so-called  "iteration-series" 
approach  already  mentioned  requires  the  numerical  evaluation  of  an  inte- 
gral over  the  range  (-i,  i).   If  the  integration  over  this  range  is  broken 

up  into  integrations  over  the  ranges  ( -«o,  ~-0>  (i,  °°)  and  (-°o,   °o) ,  and 

(2)    i 
if  the  Hankel  function  H   (271   -r->  can  be  represented  by  its  "large- 

o      rt 

»  i 

argument   approximation  for  fairly  small  values  of  -«->  then  the  radiation 
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pattern  can  be  expressed  in  terras  cf  Fresnel  integrals.   One  of  the  advan- 
tages of  this  "Fresnel  integral''  solution  is  that  the  integrals  have  asymp- 

o       o 
totic  expansions  valid  for  values  of  <p   not  too  close  to  0   or  180  .   For 

very  thin  dielectric  the  asymptotic  formulae  show  that  a  dielectric  strip 
can  be  represented  approximately  by  two  line  sources  in  free  space  in  posi- 
tions corresponding  to  the  strip  edges  (at  x  -  +  i ) .   This  will  be  referred 
to  as  the  "two-source  equivalent". 

As  mentioned  above,  the  accuracy  of  the  Fresnel  integral  solution 

(2) 

depends  on  how  closely  the  Hankel  function  H  '  (2Mx)    is  approximated  by 

its  large- argument  asymptotic  form.   The  Hankel  function  and  its  approxi- 
mation are  compared  in  Figure  5  and  the  approximate  formula  error  is  less 
than  .015  for  x  '£  .5.   Thus  the  Fresnel  integral  solution  and  the  intera- 
tion  series  solution  should  differ  by  a  negligible  amount  as  long  as  the 
dimension  i  is  at  least  a  half  wavelength. 


2.   EXPERIMENTS 

For  convenience  in  making  the  measurements,  the  experiments  were  set 
up  on  a  ground  plane  as  shown  in  Figure  3  and  the  field  measurements  were 
made  by  using  the  corner  monopole  for  transmitting  and  the  center  monopole 
for  receiving.   A  bolometer  was  used  as  the  detecting  element  and  most 
of  the  measurements  were  carried  out  at  5000  mc .   All  the  radiation  pat- 
terns show  relative  values  of  the  electric  field  intensity  E. 

Experimentally  it  was  found  that  the  slab  height  (h)  had  negligible 
effect  on  the  radiation  pattern,  provided  that  the  height  was  greater  than 
about  three  wavelengths.   Accordingly  a  height  of  four  wavelengths  was 
chosen  for  the  experiments  performed.   The  fact  that  the  antenna  pattern 
is  independent  of  h  (for  h  sufficiently  large)  is  a  partial  justification 
for  the  use  of  a  mathematical  model  with  h  equal  to  infinity.   Probably 
the  radiation  pattern  is  independent  of  h  because  (in  free  space  at  least) 
the  monopole  radiates  no  energy  in  the  z-direction.   This  means  that  the 
top  edge  of  the  dielectric  slab  would  intercept  a  relatively  small  portion 
of  the  total  radiated  energy  if  h  were  large. 

The  theory  indicates  that  the  pattern  perturbations  depend  on 
(€   -  1)  -S- .   Experimentally  this  result  was  verified  by  selecting  various 

combinations  of  frequency,  sheet  thickness  and  dielectric  constant  to  give 

t  i 

the  same  value  of  (€   -  1)  ■=-  while  the  sheet  length  -r  was  held  constant. 

The  patterns  obtained  are  not  shown  in  this  report  but  they  were  essen- 
tially identical  for  thin  sheets  (t  <  .03^-)  and  roughly  similar  even  for 
the  thickest  sheets  considered. 


8 
The  radiation  patterns  exhibit  some  interesting  characteristics,  the 

most  prominent  of  which  is  the  "end  fire"  lobe  at  <p  =  0"'  (see  Figures  6 

to  11).   It  is  always  present  and  becomes  narrower  as  the  length  y  is 

increased.   In  addition  to  the  major  lobe  some  side  lobes  appear,  their 

i 
number  and  position  depending  on  y.   Let  us  denote  a  free  space  wavelength 

by  X-  and  the  wavelength  of  the  lowest  order  TE  mode  propagated  along  an 

infinite  dielectric  sheet  by  X-  .   It  will  be  assumed  that  X  is  approxi- 

g  g 

mately  the  wavelength  of  waves  propagating  along  a  finite  dielectric  sheet 
and  generated  by  a  line  source  as  in  the  practical  problem  considered, 
For  the  case  of  thin  sheets  X  =  X  (see  Appendix  E)  and  the  radiation 

patterns  indicate  that  the  number  of  minima  in  a  180   sector  (0  ;-:.  (p  <  180  ) 

X 

increases  by  one  for  every  —  increase  in  i .   This  rule  breaks  down  for 

the  thicker  sheets  for  which  X  zk  X  (see  Figures  10  and  11  for  the  — "  and 

g  8 

— "  sheets)  but  the  rule  seems  to  hold  for  all  thicknesses  considered  if 

we  write  X  in  place  of  X.   Thus  the  rule  is  that  the  number  of  nulls 

g  X 

o  o  fir 

in  the  range  0  <[  (p  <Z   180   increases  bv  one  for  every  ~~   increase  m  the 

X  2 

length  i.      Values  of  -e*  for  various  sheet  thicknesses  are  tabulated  in 

Appendix  E  where  it  can  be  seen  that  X   is  always  somewhat  less  then  X, 

X  E 

Since  the  values  of  -j—   decrease  rapidly  for  t  >  „05292X  (•-"*)  it  is  not 

surprising  that  the  effects  of  low  propagation  velocity  show  up  in  the 
radiation  patterns  only  for  t  =  .1058X  (—" ) . 

Another  interesting  aspect  of  the  radiation  patterns  is  the  null  depth, 
a  quantity  whicn  increases  sharply  as  the  slab  thickness  and/or  dielec- 
tric constant  increases.   In  all  cases  the  null  beside  the  major  lobe 
(which  we  shall  call  the  principal  null)  is  the  deepest  while  all  the 
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15 
other  nulls  have  roughly  equal  depth.   This  immediately  suggests  a  means 
of  expressing  the  magnitude  of  perturbation  of  the  original  circular  dipole 
radiation  pattern  which  is  as  follows: 


Let    P  =  unperturbed  dipole  pattern 
(not  a  function  of  (p) 

A  P  =  maximum  pattern  perturbation  or  error 

P    =  magnitude  of  major  lobe 
max     &  J 


P  .   =  magnitude  of  principal  null 
mm 


Thus  A  P  is  defined  to  be 


P    -  P 

A  p    max    min 


In  this  report  P  is  taken  to  be  unity  so  that 


A  P  =  P    -  P 

max    mm 


In  addition,  it  is  convenient  to  express  A  P  as  a  percentage.  The  experi- 
mental values  of  A  P  are  given  in  Figures  12  to  17  where  they  are  compared 
with  the  theoretical  values  obtained  from  the  "iteration  series"  solution. 
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3 .   THEORY 

The  problem  of  the  scattering  of  electromagnetic  waves  by  a  dielec- 
tric body  is  a  very  difficult  one  to  solve  and  in  general  approximate 
solutions  must  be  sought.   Since  approximate  methods  are  especially  use- 
ful in  the  solution  of  integral  equations,  it  was  decided  to  formulate 
the  problem  considered  in  this  report  as  an  integral  equation.   The  method 
of  induced  sources  provides  a  convenient  means  of  deriving  this  equation 
and  also  gives  some  insight  into  the  manner  in  which  the  dielectric  affects 
the  free-space  radiation  pattern  of  the  metallic  antenna. 

The  problem  of  scattering  from  an  arbitrary  inbomogeneous  region  in 
space  can  be  formulated  as  follows.   Consider  the p roblem  of  determining 
the  total  electromagnetic  field  produced  by  sources  J  and  K  (electric  and 
magnetic  currents)  in  free  space  in  the  presence  of  an  inbomogeneous  region 
R.   In  free  space,  Maxwell's  equations  have  the  form 


curl   H  =  Y     E   +  J  Y=jw€ 

o  o 


curl    E   = -Z      H-K  Z      =   j    w  U 

o  o  o 


and  in  region  R  they  have  the  form 


curl  H  =  Y„  E      Y  =.  j  w  £  £ 
R         R   J     o  r 


curl  E  ^-Z„  H      Z  =  j  u  fi  |i 
R         R   J    *o  *r 


(If  polarization  or  conduction  losses  in  R  cannot  be  neglected  €   and  U 

fe        r      r 

may  be  complex.)   In  order  to  achieve  a  formulation  in  infinite,  homogene- 
ous free  space,  the  equations  for  region  R  can  be  rewritten  as 
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curl  H=Y   E  +  (Yn-  Y  )  E 
o       R     o 


curl  E  =  -Z   H  -  (Z„  -  Z  )  H 
o       R     o 


and  the  terms  containing  the  differences  in  medium  parameters  can  be  inter- 
preted as  source  terms  or  "induced  sources".   Thus  the  total  field  can 
be  constructed  by  a  superposition  of  the  fields  due  to  two  sources  in  free 
space.   The  incident  field  is  due  to  the  primary  sources  J  and  K  and  the 

scattered  field  is  due  to  the  induced  sources  J.  and  K.  where 

1      1 

J,  =  (Y  -  Y  )  E 

l      R     o 

K.  =  (Z„  -  Z  )  H 
l      R     o 

This  procedure  does  not  lead  directly  to  a  solution  of  the  problem  since 
E  and  H  are  unknown  in  region  R.   However  it  does  permit  the  formulation 

of  the  problem  as  an  integral  equation  for  the  total  field  as  shown  by 

2 

Rhodes   who  considers  the  case  of  a  lossless  dielectric  scatterer  in  which 

ji  =  1.   In  this  case  the  induced  current  is  entirely  electric  and  has 
the  form 


J.  =  j  w  €  (€   -  1)  E 
l        or 


It  is  now  possible  to  see  how  the  dielectric  scatterer  perturbs  the 
field  of  a  metallic  antenna.   If  the  induced  sources  are  distributed  over 
a  region  in  space  having  at  least  one  dimension  comparable  to  a  wavelength, 
then  they  will  contribute  a  significant  amount  of  radiated  field  which 
is  not  present  if  the  metallic  antenna  is  located  in  free  space.   In  such 
a  case  the  "total"  antenna  is  composed  of  dielectric  as  well  as  metal  parts 
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2 
The  integral  equation  can  be  specialized  (as  shown  by  Rhodes  )  for 

the  case  of  an  infinite,  cylindrical  dielectric  scatterer  of  uniform  but 
arbitrary  cross-section  when  the  incident  wave  propagates  in  a  direction 
normal  to  the  cylinder  axis  and  is  polarized  parallel  to  the  axis.   Appen- 
dix A  shows  how  this  integral  equation  can  be  modified  for  the  case  of 
an  infinite  line  source  centrally  located  inside  an  infinite  dielectric 
cylinder  having  the  cross-section  of  a  long,  thin  rectangle.   The  solu- 
tions of  this  integral  equation  are  derived  in  detail  in  the  appendices 
and  summarized  in  the  remainder  of  this  chapter. 

The  Iteration  Ser 1 e s  Solution 

The  details  of  this  method  are  presented  in  Appendix  A.   In  the  deri= 
vation  it  is  assumed  that  the  dielectric  sheet  is  thin,  that  its  relative 
dielectric  constant  is  low,  that  there  is  negligible  field  variation  across 
the  narrow  dimension  of  the  sheet,  and  that  the  field  in  the  sheet  is 
the  same  as  the  incident  field  in  free  space. 

Since  the  experimental  dipole  was  inserted  in  a  small  hole  drilled 
into  the  edge  of  the  dielectric  sheet,  it  is  necessary  to  consider  a  line 
source  which  is  not  embedded  in  the  dielectric  but  rather  located  centrally 
in  a  small  hole  of  radius  6.   Both  experimental  and  theoretical  radiation 
patterns  seemed  quite  insensitive  to  changes  in  6  provided  that  it  was 
kept  very  small  compared  to  a  wavelength.   For  most  of  the  experiments 
described  in  this  report  the  central  hole  had  a  radius  of  about  a  hundredth 
of  a  wavelength  so  the  value  6  =  „OlX.  was  chosen  for  the  calculations. 

The  derivation  gives  the  following  radiation  pattern  formula  from 
which  all  the  theoretical  results  in  Figures  6  to  17  were  calculated, 


where 
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P(0    =  \J  1    -    4KB  +  4K2  (A2  +  B2) 


i 

A  =      J  (27Tx)  cos  (27Tx  cos  <£>)  dx 
.01 


i 

."X 

B  =      N  (21Tx)  cos  (21Tx  cos  0)    dx 
J     ° 
,01 


K  =ff2(€r  -  1)  * 


The  integrals  A  and  B  are  tabulated  in  Appendix  B,  together  with  numeri- 
cal values  of  P(<jP)  . 

The  above  formula  for  the  radiation  pattern  can  be  simplified  even 
further  for  the  case  of  extremely  thin  dielectric  sheets.   For  this  case 
K  «  1  and  K2  «  K  which,  together  with  the  fact  that  A  and  B  are  of  the 
order  of  unity,  gives  tne  following  approximation  for  the  radiation  pat- 
tern. 

P(Cp)    =  1  -  2KB 

Thus,  for  extremely  thin  dielectric  sheets,  the  pattern  perturbations 

t 
are  directly  proportional  to  (C   -  1)  y.   Furthermore  P((0    is  dependent 

only  on  B,  an  integral  involving  the  Neumann  function  N  .   There  is  only 

a  secondary  contribution  from  A,  an  integral  involving  the  Bessel  function 

J  . 
o 
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The  derivation  summarized  above  (and  described  fully  in  Appendix 
A)  gives  the  field  at  large  radius  due  to  a  line  source  at  the  origin. 
An  identical  radiation  pattern  formula  can  be  derived  in  somewhat  fewer 
steps  by  applying  the  principle  of  reciprocity.   This  means  that  the  radia- 
tion pattern  is  given  by  the  magnitude  of  the  electric  field  at  the  origin 
due  to  an  incident  plane  wave  coming  from  sources  at  infinity. 

The  Fresnel  Integral  Solution 

This  method  of  solving  the  Integral  equation  is  presented  in  detail 
in  Appendix  C.   The  assumptions  involved  are  the  same  as  those  in  the 
iteration  series  solution  with  two  exceptions.   First,  in  the  Fresnel 
integral  solution  it  is  assumed  that  the  central  line  source  is  embedded 
in  the  dielectric.   This  is  done  in  order  to  take  advantage  of  a  very 
simple  integral  formula  and  is  in  contrast  with  the  iteration  approach 
which  considers  the  line  source  to  be  in  a  small  hole  of  radius  6.   Second, 

it  is  assumed  that  the  large-argument  approximation  for  the  Hankel  func- 

(2)    i 
tion  H"   (27T  y)  is  accurate  even  for  fairly  small  arguments.   Figure  5 

shows  this  assumption  to  be  valid  since  the  approximate  formula  error 

is  less  than  .015  if  i  >  ,5X. 

The  basic  idea  m  the  Fresnel  integral  approach  is  that  of  changing 

the  range  of  integration  from  (-i,  i)  to  a  combination  of  (-00,  -i),  (£,    °o) 

and  (-=oc,  00),   This  is  equivalent  to  breaking  down  the  field  scattered 

by  a  finite  dielectric  strip  (from  x  =  -i  to  x  =  i)  into  the  difference 

between  the  field  scattered  by  an  infinite  sheet  and  the  field  scattered 

by  two  semi-infinite  strips  (from  x  =  -*©  to  x  =  -i  and  from  x  =  I    to 
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x  =  oo) .   This  viewpoint  gives  a  solution  to  the  integral  equation  from 
which  the  radiation  pattern  P(<#)  can  be  obtained. 


H<P)    = 


1  -  7T(€   -  1)  y  •  — ^ 
r      X   sin  <f> 


M  +  j(l  -  N) 


where 


M  =  sin  -r 

mm 


N  =  sin  % 


1-S 


2\  / 

1 

irX 


-  c 


/2V  i 

/lfx~ 


-  c 


irx 


7TX 


+  cos  - 


+  cos  - 


2V  i 

2 


-  C 


1-S 


27  i 

2 

7TX 


2y  I 


-  C 


V      =    277(1    +  cos   <?),       "Y      =   277(1    -    cos   (?) 
1  2 

7  V       V 

Tables   of  S  and  C  are  available  and  give  S  and  C  as  functions  of  — r— and  7    • 

A.      K 

Thus  the  radiation  patterns  can  be  calculated  in  a  straightforward  manner 
without  the  aid  of  a  machine  computer. 

Results  of  the  iteration  and  Fresnel  integral  calculations  are  compared 
in  Figure  18  for  I    =  2.0X  and  t  =  .02646X.   The  pattern  magnitudes  given 
by  the  two  methods  differ  by  a  small,  almost  constant  amount  (.035  approxi- 
mately).  This  difference  is  believed  to  be  due  to  the  fact  that  the  itera- 
tion approach  considered  the  line  source  to  be  in  a  hole  of  radius  6  =  .OlX, 
whereas  the  Fresnel  integral  solution  takes  6  to  be  zero.   This  explana- 
tion was  substantially  verified  by  some  approximate  calculations  (itera- 
tion series  method)  which  indicated  that  a  pattern  difference  of  about 
.03  was  to  be  expected. 


FIGURE  18 

Comparison  of  the  Iteration  Series,  Fresnel  Integral 
and  "Two-Source  Equivalent"  Radiation  Patterns 
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Asymptotic  Behavior  of  the  Fresnel  Integral  Solution 

Appendix  D  shows  that,  when  <p   is  not  too  close  to  0   or  180  ,  the 
field  at  the  origin  P  due  to  an  incident  plane  wave  is  given  by 


t  f-j 

E(P)    =    1    +  7T(€      -   1)  v    1  ™^77>  +  -o 
r  A.        sin  (p         2 


-J27T  £  cos   <p  +JZII  ~   cos   (p. 

e  e 

-  + 


R<|) 


1  +  COS  <^ 


1  -  COS  <? 


where  R(-r)  is  the  large-argument  approximation  to  the  Hankel  function 

(2 )    £  i 

H   (271^-).   The  factors  involving  R(^r)  behave  as  if  they  were  due  to  infi- 
nite line  sources  at  x  =  -i  and  x  =  +1    with  directional  characteristics 

-I  v-1 

(1  +  cos  (p)        and  (1  -  cos  <p)        respectively  and  with  phase  factors 

i  A 

-j27f  T  cos  Cp  +J27T  t  cos  <f> 

e     *>       and  e  respectively. 

Some  further  restrictions  lead  to  a  simple  formula  for  the  radiation 

pattern  (^the  absolute  value  of  E(pV)  ,   If  (C   -  1)  -=-  is  very  small,  the 

—1 
term  (j  sin  (p)        gives  a  negligible  contribution  to  the  radiation  pattern 

and  only  the  line  source  terms  remain.   Furthermore,  if  i  is  some  integral 

multiple  of  a  half  wavelength  (i  =  n  —  as  in  the  patterns  in  Figures  6 

to  11),  the  radiation  pattern  formula  becomes 


t         f-1) 
P((p)  =  l  +  (€     -  i)  --   .  A_LL 

r  ^  f— 


cos(nXT   cos  (0    -   cos   <p  sin(n7T   cos  (p) 
sin2    (p 


This   pattern   is   plotted   in   Figure   18   for  i    =   2. OX  (n  =   4)    and   t   =    .02646X. 
Numerical   values   are   as    follows: 
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(p  P((?) 


10° 

1.795 

20° 

1.241 

30° 

1.061 

46° 

.952 

62° 

1.029 

76° 

.978 

90° 

1.021 

In  the  above  expression  for  P(.(f) ,    the  line  source  directivity  factors 

-1 

(1  +  cos  <P)    "    have  been  retained  and  the  formula  holds  for  values  of  <p 

o 
as  small  as  20  . 

The  formula  for  P(<p)    can  be  simplified  even  further  by  considering 
the  line  sources  to  be  non-directive.   This  means  that  cos  (f *  =   0  and 
sin  <p  =  1 .   Thus  <p   is  restricted  to  be  fairly  close  to  90  .   The  radia- 
tion pattern  formula  becomes 


t    (-l)n 
P((p)  =  1  +  (£   -  1)  T  •  v  '      cos(n1T  cos  (p) 

/n 


This  pattern  is  also  shown  in  Figure  18  for  i  =  2.0X  (n  =  4)  and 

t  =  .02646X.   Numerical  values  are  as  follows  (they  are  taken  at  pattern 

maxima  and  minima): 

(p  P((?) 


0° 

1.021 

43° 

.980 

62° 

1.021 

76° 

.980 

90° 

1.021 
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In  the  above  expression  for  P((p)    the  directivity  of  the  individual  line 
sources  has  been  neglected.   Despite  this  severe  restriction,  the  formula 

holds  fairly  well  for  values  of  cp  as  small  as  40  .   In  particular,  the 

o 
angles  at  which  minima  and  maxima  occur  agree  to  within  +  3  with  both 

the  experimental  results  and  the  iteration  series  calculations  for  very 
thin  dielectric  sheets. 

The  preceding  discussion  of  the  "two-source  equivalent"  has  been 
presented  because  it  seems  to  bring  out  two  factors  which  help  in  under- 
standing the  phenomena  described  in  this  report.   First,  as  far  as  radia- 
tion pattern  calculations  are  concerned,  it  is  more  important  to  know  the 
sheet  length  i  than  the  exact  field  distribution  along  the  sheet.   Second, 
one  can  think  about  the  problem  in  terms  of  outward-traveling  cylindrical 
waves  which  are  scattered  at  the  sheet  edges  and  nowhere  else. 
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4 .   COMMENTS 

There  is  close  agreement  between  the  "iteration  series"  theory  and 
the  experiments  from  t  =  0  to  t  =  .03X  (for  €  =  2.55),  within  which  range 
the  pattern  perturbations  vary  almost  linearly  with  respect  to  sheet  thick- 
ness.  As  t  increases  beyond  ,03X  the  theory  becomes  useless  in  that  it 
does  not  predict  the  observed  radiation  patterns,  especially  the  deep  nulls 
Such  a  breakdown  of  the  theory  is  not  surprising  since  most  of  the  assump- 
tions involve  a  very  small  thickness  and  a  low  dielectric  constant.   Im- 
proved results  probably  could  be  obtained  by  computing  one  or  more  addi- 
tional terms  in  the  iteration  series  solution  for  the  field  in  the  dielec- 
tric sheet. 

The  results  reported  for  €  s  2.55  can  be  extended  easily  to  other 
values  by  use  of  the  fact  that  the  pattern  perturbations  do  not  depend 

on  t  alone  but  rather  on  (€   -  1)  ■=-.   This  dependence  is  complicated  for 

X. 

moderately  thick  dielectric  sheets  (of  the  order  of  —  thick)  but  is  prob- 
ably of  the  form 


>  X 


r      a.  j  i 


I   + 


(€   -  1)  ± 
r      X 


2 

I   + 
2 


I 
where  1^  and  I   are  factors  depending  on  cp  and  y-   This  form  is  suggested 

by  the  infinite  series  solution  which  can  be  derived  by  applying  the  itera- 

i 
tion  technique  to  the  integral  equation.   For  constant  length  -^,    dielectric 

sheets  having  the  same  value  of  (€   -  1)  -^  should  produce  nearly  the  same 

radiation  pattern.   Experimentally  this  was  found  to  be  true  in  the  linear 

range  where  (€   -  1)  r  is  very  small  and  approximately  true  even  for  the 

thickest  dielectric  considered  in  this  report. 
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For  thin  dielectric  sheets  it  is  convenient  to  express  the  radiation 
pattern  perturbation  or  error  as 


P    -  P 

A  _    max    mm 
A  P  =  — — — - 


where  P  is  the  magnitude  of  the  unperturbed,  circular  dipole  radiation 
pattern.   The  results  obtained  for  pattern  error  are  higher  than  one 
might  expect;  that  is  a  A  P  of  14%  occurs  when  a  dipole  is  laid  on  a 
square  polystyrene  sheet  six  wavelengths  on  a  side  and  only  a  hundredth 
of  a  wavelength  thick „ 

The  theory  described  in  this  report  is  helpful  in  obtaining  inter- 
pretations of  the  observed  radiation  patterns.   For  example,  the  view- 
point of  induced  sources  results  in  a  useful  integral  equation  and  per- 
mits the  dielectric  to  be  regarded  as  a  distribution  of  electric  currents 
In  free  space.   This  concept  leads  to  the  conclusion  that  a  dielectric 
object  in  the  near  field  of  an  antenna  will  alter  the  radiation  pattern 
perceptibly  if  its  size  is  comparable  to  a  wavelength  in  at  least  one 
dimension „ 

Another  viewpoint  is  that  of  outward- traveling  cylindrical  waves 
which  are  scattered  only  at  the  sheet  edges.  This  concept  suggests  that 
the  dielectric  sheet  may  be  represented  approximately  by  line  sources 
in  free  space  in  positions  corresponding  to  the  sheet  edges.  (Such  a 
representation  is  well  known  for  the  problem  of  scattering  by  a  metallic 
strip).  The  asymptotic  approximation  to  the  Fresnel  integral  formula 
justifies  such  a  representation  when  the  sheet  is  very  thin  and  when 
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the  radiation  field  is  measured  in  a  sector  roughly  perpendicular  to 
the  broad  face  of  the  sheet  (<p   is  not  close  to  0   or  180  )  .   A  more 
accurate  representation  results  if  the  individual  line  sources  are  con- 
sidered to  have  non-circular  radiation  patterns. 

All  the  calculations  in  this  report  are  for  a  dipole  source  located 
half  way  between  the  edges  of  the  dielectric  sheet.   This  approach  sim- 
plifies the  theoretical  derivations  but  the  formulae  for  the  radiation 
pattern  of  a  dipole  at  any  point  on  the  sheet  could  be  derived  in  a 
similar  way.   In  particular,  the  Fresnel  integral  radiation  pattern 
formula  can  be  altered  easily  to  fit  the  non-symmetrical  problem.   If 

the  dielectric  strip  extends  from  x  =  -i   to  x  =  +1      and  if  the  source 

1  2 

remains   at   x  =   0,    then   the   radiation  pattern  is 


P(<£)    =      l    -  7T<€      -   1)  1    .  -J— 
r  X        sin  (p 


M   +   j(l    -   N) 


where 


sin* 

f 

|2V  i      \ 

f 

/2\  i      \ 

/ 

/2\  i      \ 

( 

/2Y  4     \ 

M   = 

K 

V  **     ) 

-c( 

V  irx     ) 

+   cos  - 

1                    • 

s(v 

'tPJ 

-cL 

hrPjj 

sir^ 

/ 

<  lz\  t   > 

i               / 

'  [2Y  t     \ 

+    COSr 

r  jzy  t    N 

\ 

/  /2Y  i 

N  = 

i-( 

/     i  i 

J     TTX       j 

)-c( 

1      11 

1-S 

W     TTX        , 

)-<= 

Wie 

Y      =    2JT(1    +   cos   <£) 

1 


Y     =    2JT(1    -    cos   <£) 

2 
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5.   CONCLUSIONS 

The  experimental  radiation  patterns  have  been  shown  for  the  case 
of  a  dipole  embedded  inside  a  dielectric  sheet  (or  "printed"  on  its  sur- 
face) and  these  results  have  been  found  to  agree  with  the  theoretical 
calculations  for  thin  sheets.   The  range  of  thickness  over  which  theory 
and  experiment  agree  is  a  useful  range  from  the  point  of  view  of  practi- 
cal "printed  circuit"  antennas.   As  a  consequence  the  results  reported 
permit  tne  antenna  designer  to  estimate  the  maximum  radiation  pattern 
error  likelv  to  be  encountered  wnen  mounting  single  dipoles  or  combina- 
tions of  them  on  dielectric  sheets. 

The  theory  presented  shows  that  a  very  thin  dielectric  strip  can 
be  regarded  as  being  approximately  equivalent  to  two  line  sources  in 
free  space  in  positions  corresponding  to  the  strip  edges.   This  result 
i?.    quite  general  and  could  be  applied  to  many  problems  involving  the 
scattering  of  electromagnetic  waves  by  very  thin  dielectric  sheets „ 
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APPENDIX  A 
The  Iteration  Series  Solution  of  the  Integral  Equation 

For  the  case  of  a  dielectric  scatterer,  the  integral  equation  for 

2 
the  two-dimensional  scattering  problem  as  given  by  Rhodes   is 


where 


Ck  2  -  k  2) 

E  (P)  =  E  X(P)  -  j  -i ° 

z       z  4 


E  (Q)  H  (2)(k  p)  da 

zoo 


A 


P  =  (x,  y)  or  (r,  <p)    and  is  the  observation  point 
Q  =  (x'>  y')  or  (r'  <p  )  and  is  a  general  point  inside 
the  dielectric 

/N2 


P  ~  /(x  -  x')2  +  (y  -  y7) 


r  =   /  2  ,   2 
vx  +  y 

27T 
k  =  -ST-  =  the  propagation  constant  in  free  space 

0  A. 

27T   / 

k  =  -r-  /€  =  the  propagation  constant  in  the  dielectric 

1  A,  \  r 

€   =  the  relative  dielectric  constant  of  the  scatterer 
r 

X  =  the  free  space  wavelength 

A  =  the  cross-sectional  area  of  the  scatterer 

\ (P)  =  the  total  electric  field  at  P„   (E  =  E  =  0) 

E  (P)  =  the  incident  field  at  P  due  to  a  two-dimensional 
z 

source  in  free  space 
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The  above  general  equation  must  be  specialized  for  the  case  of  an 
infinite  line  source  located  centrally  inside  an  infinite  dielectric 
strip  as  shown  in  Figure  4.   As  discussed  in  Chapter  4,  it  is  necessary 
to  consider  the  line  source  as  being  in  a  hole  of  radius  6  since  the 
experimental  dipole  was  inserted  in  such  a  hole  drilled  into  the  edge 
of  the  dielectric  sheet,   It  will  become  evident  that  such  a  considera- 
tion avoids  integrating  over  the  Hankel  function  singularity  at  the  origin. 
In  the  derivation  to  follow  it  will  be  understood  (if  not  explicitly 
stated)  that  any  integration  over  the  cross-sectional  area  of  the  dielec- 
tric sheet  excludes  a  circular  region  of  radius  6  around  the  line  source 
at  the  origin. 

If  it  is  assumed  that  the  dielectric  sheet  is  very  thin  and  that 
there  is  no  field  variation  across  it,  then 


E  (Q)  H<2)(k JO)  da  =  t     E  (Q)  H(2)(k  p)  dx' 

z     o    o  J    z     o    o 

-1 


where  t  is  the  sheet  thickness  and  2i  is  its  length.  Using  this  relation 
and  writing  E  in  place  of  E  ,  we  obtain  an  approximate  expression  for  the 
total  field  which  is 


(k   -  k   )     r 

E(x,  y)  =  C  H(2)(k  r)  -  j  ~l—~-^~-   t     E(x',  o)  H(2>(k  p)  dx/ 

o     o  4         J       '  o     o 

-i 

(?) 

where  C  H  "  (k  r)  is  the  field  of  an  infinite  line  source  located  at 
o    o 

the  origin,  C  being  a  constant  proportional  to  current  strength. 
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Letting  y  =  0  in  the  above  expression  gives  an  approximate  one- 
dimensional  integral  equation  for  the  field  in  the  dielectric  sheet. 

(    \  (k-k)r  (    . 

E(x,  o)  =  C  r2\k  x)  -  j  — i A      °        t      E(x',  o)  H^2'  (k  |x  -  x'l)  dx/ 

7  o     o  4  o      o  ■       ■ 

-I 

If  this  equation  is  solved  then  the  total  field  can  be  obtained  from  the 
above  expression  for  E(x,  y) .   For  the  solution  of  this  integral  equation 
the  iteration  technique  is  to  be  employed.   It  consists  of  making  an 
approximation  to  E(x/,  o),  calculating  E(x,  o),  and  reinserting  this 
function  under  the  integral  sign  in  place  of  E(xx,  o) .   By  repeating 
this  process,  successively  more  accurate  values  of  E(x,  o)  can  be  obtained, 
In  other  words,  a  series  expression  for  E(x,  o)  can  be  derived  and  this 
series  is  known  to  converge  if  the  coefficient  preceding  the  integral 

is  small  enough.   For  these  calculations  it  will  be  assumed  that 

/  (2) 

E(x  ,  o)  =  0  which  gives  E(x,  o)  =  C  H   (k  x) .   The  iteration  procedure 

00 

will  be  carried  no  further  than  this.   That  is,  we  shall  consider  only 
the  first-order  term  in  the  iteration  series.   This  is  simply  a  state- 
ment of  the  assumption  that  the  field  in  the  dielectric  sheet  is  equal 
to  the  incident  field  in  free  space  represented  by  the  zero-order  Hankel 
function.   This  value  of  E(x",  o)  will  be  inserted  after  some  simplifi- 
cation of  the  formula  for  E(x,  y) . 

At  this  point  the  notation  can  be  made  somewhat  more  compact  by 

expressing  all  lengths  in  terms  of  free  space  wavelengths,  using  bars 

—    r  —    x  — 

(— )  to  indicate  the  new  coordinates  so  that  r  =  t,  x  =  v  and  k  r  =  27Tr. 

V      \      o 

1    2      o     7T2 
Thus  it  can  be  seen  that  — (k   -  k   )  =  — (€   -  1)  and  that  the  integral 

4   1      o     X2   r 
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expression  for  E(x,  y)  becomes 


2)/™,^    ^2tc  ,\  -     J   w~/  „n  tr(2)/0(Tr^N  ^-/ 


E(x,  y)  =  C  H^'(27Tr)  -  jTT *(€^  -  1)  t  J   E(x  ,  o)  IT"(27rp)  dx 


-i 

Setting 

p  =  -jir2(€  -  i)  t 

r 
and 


I  =  J       E(x',  o)  H(2)(27Tp)  dx^ 


we  can  put  the  integral  expression  for  E(x,  y)  in  the  form 

E(x,  y)  =  C  H(2)(2pr)  +  Pi 

The  integral  I  can  be  transformed  to  an  integral  over  the  range 
0  <  x  _<  i  by  the  following  method.   We  have 


I  =  J        ECx',  o)  H^2)(27Tp)  dx' 


-i 

Recalling  that  the  line  source  is  considered  to  be  in  a  hole  of  radius  8 
(or  6) ;  we  have 

-"5  1 

I  =J    E(x',  o)  H^2)(27Tp)  dx'  +J   ECx',  o)  H*2)(27Tp)  dx' 
-7  5 
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From  the  symmetry  of  the  problem  it  is  apparent  that  E(x/J  o)  =  E(-x  ,  o) 
Then,  writing  -x   for  x  ,  we  have 


I 


x.  x. 

=  J   E(-x',  o)  H(2)[27Tp(-x/)]  dx'  +J   E(x',  o)  H(2  *[  27Tp(x')  ]  dx' 


=  J   E(x',  o)   J  H^2)[27Tp(-x^)]  +  H^2)[27rp(xO]  I   dx' 

Now  a  major  simplification  follows  from  the  observation  that  p   and  r  are 
very  large  for  the  case  of  radiation  field  calculations.   Thus  the  follow- 
ing large-argument  approximations  for  the  Hankel  functions  can  be  made. 


"j(k  r  -  -) 


H<2>(kr)=  /  »    e  ^  o    ^ 


o     o      \/  7Tk  r 


o    o^     ,MTkp 


If  r  and  p  are  large,  we  can  set  p  =  r  in  the  amplitude  factor  and 
P  =  r  -  x  cos  (p   in  the  phase  factor  so  that 


J. 
-j(k  r  -  — )     jk  x  cos  cp 


H   (k  P)  =   /— e  •  e 

o    o      ,/  TTk  r 


Using  the  barred  notation,  the  Hankel  function  approximations  are  given 
by 
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(2)        -  1  -J(2*r-*) 

o 


*f 


H<2>(21rp)   .  _l_     e"J<arr   "  4}    .    eJ2>r*  cos  ? 

o 


TT  v/  r 


Introducing   the   above   largQ-argument   approximations,    we  have 


fl  t  -j(27Tr   -  ^)      r      .-,  (f)  -j27Tx'cos   <jp   , 

I   =  E(x   ,    o)   — —      e  \  e  +      e  f      dx 


7T  /7 


i*  —         TT 

r    -  ~J(27rr  -  v 

E(x',  o)  — —      e  •    2   cos( 2^x^03   (?)    dx/ 

"3  7T  /  r 


i 
=  J       E(x',     o)    N(r,    x',    (?)    dx' 
6 


where 


N(r,    x',    (p)   =   R(r)    X(x ',    <j0) 


-j(2frr   -  ^) 
R(r)    =  — —      e 


7T  >/  r 


X(x',    <p)    =    2   cos(27Tx'cos   (?) 
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The  problem  can  now  be  stated  entirely  in  a  polar  coordinate  system  in 
which  E(x,  y)  becomes  E(r,  (p)    and  E(x/,  y/)    becomes  E(r/,  <p')  .      In  simpli- 
fied form  the  integral  expression  for  the  radiation  field  is 

/ 

E(r,  (p)    =  C  R(r)  +  (3  J   E(r',  o)  N(r,  r't    (p)    dr' 

6 


The  above  relation  will  permit  a  calculation  of  E(r,  cp)  provided 
that  E(ry,  o)  is  known.  The  field  in  the  dielectric  sheet  is  unknown 
for  all  cases  except  the  case  of  an  infinitely  thin  sheet  with  a  low  dielectric 

constant.   For  that  case  the  field  in  the  sheet  will  be  approximately  equal 

( o)  / 

to  the  incident  field  in  free  space.   That  is  we  can  set  E(r  ,  o)  =  CH   '  (27Tr  ), 

o 

an  approximation  already  shown  to  be  equivalent  to  the  use  of  the  first  order 
approximate  solution  of  the  integral  equation  for  the  field  inside  the 
dielectric  sheet.   Thus  we  have 


E(r,  (p)    =  C  R(r)  +  (3  R(r)  C 


"5 


H(2)(27Tr/')  X(r'  (p)    dr' 
o  ' 


C  R(r) 


1+6    H(2)(27Tr/)  X(r'  <p)    dr' 
8 


and  the  radiation  pattern  P(<p)    is  given  by 


P(<?)  = 


1  +  Pj   H^2)(27Tr^)  X(rx,  (p)    dr' 
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-/ . 


The  variable  r  is  just  a  variable  of  integration  so  the  notation  can  be 


-/ 


simplified  by  writing  x  in  place  of  r  .   The  pattern  expression  becomes 


P(<?)   = 


1    +  (3  J       H(2)(27Tx)    X(x,    CO)    dx 


1+2(3 


J    (27Tx)    cos(27Tx  cos   CO)    dx 
J  o 


"   JJ       No 


(27Tx)    cos(27Tx   cos   CO)    dx 


Setting 


A  =  j       J    (27Tx)    cos(2flx  cos   CO)    dx 


we  have 


B  =  N    (27Tx)    cos(2flx   cos   CO)    dx 

J  o 

"5 


P(<0>    =    I  1    +   2p(A   -   jB) 


Setting  (3  =   -   jK  and  K  =  7T  (€      -   1)    t,    we  have 


P(<^)  = 


1  -  2KB  -  J2KA   =   j/(l  -  2KB)2  +  (2KA)2  =   /l  -  4KB  +  4K2  (A2  +  B2 ) 


This  is  the  required  approximate  radiation  pattern  formula.   The  value 

of  6  was  set  at  .01  since  in  the  experimental  work  the  monopole  was  inserted 

in  a  hole  having  a  radius  of  about  a  hundredth  of  a  wavelength.   The  above 
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results  are  summarized  in  Chapter  3  and  numerical  calculations  are  tabu- 
lated in  Appendix  B, 
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APPENDIX  B 

Numerical  Results 

The  numerical  results  of  the  radiation  pattern  calculations  are  con- 
tained in  the  following  pages.   The  calculations  given  are  for  the  itera- 
tion series  solution  of  the  integral  equation  and  the  polar  graphs  of 
these  values  are  shown  in  Figures  6  to  11. 

Also  included  are  numerical  values  for  the  integrals 


i 

X 


A  = 


J 
.01 


J    (27Tx)    cos(27Tx   cos   cp)    dx 


i 


B   = 


J 
.01 


N   (27Tx)    cos(27Tx   cos   (0)    dx 
o 
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Radiation  Patterns  -  Table  1 


t  = 

t   =   .5X 

.01323X 

t  =  i.ox 

€r 

=   2.55 

It    =   2.5X 

f 

i    =   1.5X 

I    =   2. OX 

i    =   3.0X 

0 

1.057 

1.083 

1.103 

1.121 

1.136 

1.150 

2 

1.056 

1.083 

1.103 

1.121 

1.136 

1.150 

4 

1.056 

1.083 

1.103 

1.120 

1.135 

1.148 

6 

1.056 

1.082 

1.102 

1.118 

1.132 

1.145 

8 

1.056 

1.081 

1.100 

1.116 

1.129 

1.141 

10 

1.055 

1.080 

1.098 

1.113 

1.125 

1.135 

12 

1.055 

1.079 

1.096 

1.109 

1.120 

1.128 

14 

1.054 

1.077 

1.093 

1.104 

1.113 

1.119 

16 

1.054 

1.075 

1.089 

1.098 

1.105 

1.108 

18 

1.053 

1.073 

1.085 

1.092 

1.095 

1.095 

20 

1.052 

1.070 

1.080 

1.084 

1.084 

1.080 

22 

1.051 

1.067 

1.074 

1.075 

1.071 

1.063 

24 

1.050 

1.064 

1.068 

1.065 

1.057 

1.045 

26 

1.048 

1.060 

1.060 

1.054 

1.042 

1.027 

29 

1.047 

1.056 

1.053 

1.042 

1.026 

1.008 

30 

1.045 

1.051 

1.044 

1.030 

1.011 

.991 

32 

1.043 

1.046 

1.036 

1.017 

.996 

.977 

34 

1.042 

1.041 

1.026 

1.005 

.983 

.966 

36 

1.040 

1.035 

1.017 

.994 

.973 

.960 

38 

1.037 

1.030 

1.008 

.984 

.966 

.959 

40 

1.035 

1.024 

.999 

.975 

.962 

.963 

42 

1.033 

1.017 

.991 

.969 

.963 

.971 

44 

1.030 

1.011 

.984 

.966 

.966 

.980 

46 

1.027 

1.005 

.978 

.965 

.973 

.990 

48 

1.024 

.999 

.973 

.967 

.981 

.996 

50 

1.022 

.994 

.970 

.971 

.988 

.999 

52 

1.019 

.989 

.969 

.976 

.995 

.998 

54 

1.015 

.984 

.969 

.983 

.998 

.993 

56 

1.012 

.980 

.971 

.989 

.999 

.986 

58 

1.009 

.977 

.974 

.994 

.996 

.979 

60 

1.006 

.975 

.978 

.998 

.991 

.974 

62 

1.003 

.973 

.983 

.999 

.985 

.972 

64 

1.000 

.973 

.988 

.998 

.979 

.975 

66 

.997 

.973 

.992 

.995 

.975 

.980 

68 

.994 

.974 

.996 

.991 

.974 

.986 

70 

.991 

.976 

.998 

.986 

.975 

.991 

72 

.988 

.979 

.999 

.981 

.979 

.992 

74 

.986 

.981 

.998 

.977 

.985 

.989 

76 

.983 

.985 

.997 

.976 

.989 

.984 

78 

.981 

.988 

.994 

.976 

.992 

.978 

80 

.980 

.991 

.990 

.978 

.993 

.975 

82 

A   L 

.978 

.994 

.986 

.982 

.990 

.975 

84 

.977 

.997 

.982 

.986 

.986 

.979 

86 

.976 

.999 

.979 

.991 

.981 

.984 

88 

.976 

1.000 

.977 

.993 

.977 

.988 

90 

.976 

1.000 

.977 

.994 

.976 

.990 

Radiation  Patterns  -  Table  2 
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t  =  . 

i    =    .5X 

02646X 

I    =    1.0X 

6      = 

r 

:    2.55 

l    =   2.5X 

f 

i    =    1.5X 

i    =    2.0\ 

t  =  3.0X 

0 

1.116 

1.172 

1.215 

1.252 

1.285 

1.315 

2 

1.116 

1.171 

1.215 

1.251 

1.284 

1.314 

4 

1.115 

1.171 

1.213 

1.250 

1.282 

1.311 

6 

1.115 

1.170 

1.212 

1.247 

1.278 

1.306 

8 

1.114 

1.168 

1.209 

1.243 

1.273 

1.299 

10 

1.114 

1.166 

1.205 

1.237 

1.265 

1.289 

12 

1.113 

1.163 

1.201 

1.230 

1.255 

1.276 

14 

1.111 

1.160 

1.195 

1.222 

1.243 

1.259 

16 

1.110 

1.157 

1.188 

1.211 

1.228 

1.239 

18 

1.109 

1.152 

1.180 

1.199 

1.210 

1.215 

20 

1.107 

1.147 

1.171 

1.184 

1.189 

1.187 

22 

1.105 

1.142 

1.160 

1.167 

1.165 

1.155 

24 

1.102 

1.135 

1.148 

1.148 

1.138 

1.120 

26 

1.100 

1.128 

1.135 

1.127 

1.109 

1.083 

28 

1.097 

1.120 

1.120 

1.104 

1.078 

1.045 

30 

1.094 

1.111 

1.104 

1.080 

1.047 

1.009 

32 

1.091 

1.102 

1.087 

1.055 

1.016 

.977 

34 

1.087 

1.092 

1.069 

1.031 

.988 

.952 

36 

1.083 

1.081 

1.050 

1.007 

.965 

.936 

38 

1.079 

1.070 

1.032 

.985 

.947 

.930 

40 

1.074 

1.058 

1.014 

.967 

.937 

.935 

42 

1.070 

1.046 

.997 

.953 

.935 

.948 

44 

1.065 

1.034 

.982 

.943 

.940 

.965 

46 

1.059 

1.021 

.968 

.940 

.951 

.983 

48 

1.054 

1.010 

.958 

.941 

.965 

.996 

50 

1.048 

.998 

.950 

.947 

.980 

1.003 

52 

1.042 

.988 

.946 

.957 

.992 

1.001 

54 

1.036 

.978 

.945 

.969 

1.000 

.992 

56 

1.030 

.970 

.948 

.981 

1.002 

.978 

58 

1.024 

.963 

.953 

.991 

.997 

.964 

60 

•   1.018 

.958 

.961 

.998 

.988 

.953 

62 

1.012 

.954 

.969 

1.001 

.976 

.950 

64 

1.005 

.952 

.978 

.999 

.965 

.954 

66 

.999 

.952 

.987 

.994 

.956 

.964 

68 

.994 

.954 

.99* 

.986 

.953 

.975 

70 

.988 

.957 

.998 

.976 

.955 

.984 

72 

.983 

.961 

1.000 

.967 

.962 

.986 

74 

.978 

.967 

.999 

.960 

.972 

.981 

76 

.974 

.973 

.996 

.956 

.981 

.972 

78 

.970 

.979 

.990 

.957 

.987 

.961 

80 

.967 

.985 

.983 

.961 

.987 

.955 

82 

.964 

.991 

.976 

.968 

.983 

.955 

84 

.961 

.996 

.969 

.976 

.974 

.961 

86 

.960 

1.000 

.963 

.984 

.965 

.971 

88 

.959 

1.002 

.960 

.989 

.959 

.979 

90 

.959 

1.003 

.958 

.991 

.956 

.983 

F 
t   = 
i    =    .5\ 

iadiation   Patterns    -   Table   3 

t    =   2.5X 

.05292X 
(    =   1.0X 

€ 

r 

=   2.55 

f 

1    =    1.5X 

ft    =   2.  OX 

t    =   3. 

0 

1.24 

1.361 

1.455 

1.535 

1.606 

1.671 

2 

1.241 

1.360 

1.454 

1.534 

1.605 

1.669 

4 

1.240 

1.359 

1.452 

1.531 

1.601 

1.665 

6 

1.240 

1.357 

1.449 

1.526 

1.595 

1.657 

8 

1.238 

1.355 

1.444 

1.520 

1.586 

1.645 

10 

1.237 

1.351 

1.438 

1.510 

1.573 

1.629 

12 

1.235 

1.347 

1.430 

1.499 

1.557 

1.607 

14 

1.233 

1.341 

1.421 

1.484 

1.536 

1.580 

16 

1.231 

1.335 

1.409 

1.466 

1.511 

1.546 

18 

1.228 

1.327 

1.395 

1.444 

1.480 

1.505 

20 

1.224 

1.318 

1.379 

1.419 

1.444 

1.456 

22 

1.220 

1.308 

1.360 

1.389 

1.402 

1.400 

24 

1.216 

1.296 

1.338 

1.356 

1.354 

1.336 

26 

1.211 

1.283 

1.314 

1.318 

1.300 

1.266 

28 

1.206 

1.269 

1.287 

1.276 

1.243 

1.192 

30 

1.201 

1.253 

1.258 

1.231 

1.182 

1.118 

32 

1.194 

1.236 

1.227 

1.184 

1.120 

1.046 

34 

1.188 

1.217 

1.193 

1.136 

1.060 

.984 

36 

1.180 

1.197 

1.158 

1.087 

1.005 

.936 

38 

1.172 

1.176 

1.122 

1.041 

.960 

.909 

40 

1.164 

1.154 

1.087 

.999 

.927 

.904 

42 

1.155 

1.132 

1.052 

.963 

.910 

.920 

44 

1.146 

1.108 

1.019 

.935 

.909 

.949 

46 

1.136 

1.085 

.989 

.919 

.923 

.981 

48 

1.126 

1.062 

.964 

.913 

.947 

1.008 

50 

1.115 

1.040 

.943 

.919 

.974 

1.022 

52 

1.104 

1.018 

.929 

.933 

.998 

1.022 

54 

1.093 

.998 

.922 

.952 

1.014 

1.007 

56 

1.081 

.980 

.921 

.974 

1.019 

.983 

58 

1.070 

.965 

.927 

.993 

1.012 

.956 

60 

1.058 

.952 

.938 

1.007 

.997 

.933 

62 

1.046 

.942 

.952 

1.013 

.976 

.922 

64 

1.035 

.936 

.968 

1.012 

.955 

.925 

66 

1.024 

.932 

.983 

1.003 

.937 

.941 

68 

1.013 

.932 

.996 

.989 

.928 

.961 

70 

1.003 

.936 

1.005 

.973 

.929 

.978 

72 

.993 

.942 

1.010 

.956 

.939 

.983 

74 

.984 

.950 

1.009 

.943 

.956 

.975 

76 

.976 

.960 

1.003 

.935 

.972 

.959 

78 

.969 

.971 

.994 

.934 

.983 

.941 

80 

.962 

.981 

.982 

.939 

.984 

.929 

82 

.957 

.992 

.969 

.950 

.977 

.927 

84 

.953 

1.000 

.957 

.964 

.963 

.937 

86 

.950 

1.007 

.948 

.977 

.947 

.953 

88 

.948 

1.011 

.941 

.987 

.936 

.968 

90 

.947 

l.on 

.919 

.991 

.9"51 

.974 

46 


Integrals   -   Table  4 
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l  = 

.5\ 

1  - 

=  i.ox 

<f 

A 

B 

A 

B 

0 

.14039 

-.13587 

.20639 

-.19708 

2 

.14048 

-.13575 

.20666 

-.19678 

4 

.14076 

-.13541 

.20748 

-.19585 

6 

.14121 

-.13482 

.20882 

-.19429 

8 

.14184 

-.13*01 

. 21068 

-.19208 

10 

.14265 

-.13294 

.21302 

-.18919 

12 

.14362 

-.13163 

. 21582 

-.18557 

14 

.14476 

-.13005 

.21904 

-.18120 

16 

.14606 

-.12820 

.22262 

-.17602 

18 

.14750 

-.12607 

.22650 

-.16999 

20 

.14909 

-.12364 

.23063 

-.16307 

22 

.15080 

-.12091 

.23493 

-.15523 

24 

.15264 

-.11786 

.23931 

-.14642 

26 

.15458 

-.11447 

. 24367 

-.13665 

28 

.15663 

-.11075 

. 24793 

-.12591 

30 

.15875 

-.10668 

. 25197 

-.11422 

32 

.16095 

-.10225 

.25567 

-.10163 

34 

.16320 

-.09745 

.25893 

-.08821 

36 

.16549 

-.09230 

. 26162 

-.07408 

38 

.16780 

-.08677 

.26365 

-.05937 

40 

.17013 

-.08090 

.26489 

-.04427 

42 

.17246 

-.07467 

. 26525 

-.02897 

44 

.17476 

-.06811 

.26464 

-.01371 

46 

.17704 

-.06124 

.26300 

.00122 

48 

.17926 

-.05408 

.26027 

.01556 

50 

. 18143 

-.04667 

.25644 

.02900 

52 

.18353 

-.03904 

. 25149 

.04124 

54 

.18554 

-.03123 

.24545 

.05200 

56 

.18746 

-.02329 

.23839 

.06102 

58 

.18929 

-.01529 

.23037 

.06808 

60 

.19101 

-.00727 

.22152 

.07302 

62 

.19261 

.00070 

. 21196 

.07574 

64 

.19410 

.00855 

.20186 

.07622 

66 

.19547 

.01620 

.191*0 

.07453 

68 

.19673 

.02360 

.18076 

.07081 

70 

. 19786 

.03066 

.17015 

.06532 

72 

.19888 

.03732 

.15977 

.05837 

74 

.19977 

.04350 

.14983 

.05034 

76 

.20056 

.04913 

.14052 

.04170 

78 

.20123 

.05415 

.13203 

.03289 

80 

.20179 

.05851 

.12454 

.02441 

82 

.20224 

.06214 

.11820 

.01673 

84 

.20259 

.06502 

.11312 

.01026 

86 

.20284 

.06710 

.10942 

.00537 

88 

.20299 

.06836 

.10717 

.00232 

90 

.20304 

.06878 

.10642 

.00129 

Integrals   -   Table  5 
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t 

0 
2 
4 
6 
8 
10 

12 

16 
18 
20 

22 
24 
26 
28 
30 

32 
34 
36 
38 
40 

42 

44 
46 
48 
50 

52 
54 
56 
58 
60 

62 
64 
66 
68 
70 

72 
74 
76 
78 
80 

82 
84 
86 
88 
90 


1  -- 

A 

=  1.5\ 

B 

1  = 

A 

2.  ox 

B 

.25570 
.25621 
.25773 
.26022 
.26364 
.26792 

-.24359 
- . 24304 

-.24137 
-.23855 
- . 23451 
-.22920 

.29663 

.29742 
.29977 
.30362 

.30885 
.31532 

-.28261 

-.28176 
-.27921 
-.27487 
-.26865 
-.26038 

.27295 
.27860 
. 28474 
.29117 
.29769 

-.22250 
-.21434 
-.20460 
-.19322 
-.18011 

.32279 
.33101 
.33962 
. 34824 
.35639 

-.24991 
-.23707 
-.22169 
-.20367 
-.18296 

.30407 
.31002 
.31527 
.31951 
.32245 

-.16524 
-.14863 
-.13034 
-.11051 
-.08938 

.36355 
. 36916 
.37263 
.37340 
.37092 

-.15961 
-.13382 
-.10594 

-.07653 
-.04632 

.32378 
.32321 
.32053 
.31553 
.30813 

-.06727 

-.04460 

-.02186 

.00032 

.02133 

.36479 
.35471 
.34061 
. 32264 
.30125 

-.01625 
.01257 
.03894 
.06164 
.07951 

.29829 
. 28612 
.27182 
. 25570 
. 23822 

.04046 
.05704 
.07045 
.08019 
.08591 

.27715 
.25133 
.22502 
.19958 
.17639 

.09163 
.09740 
.09674 
.09007 
.07842 

.21990 
.20135 
.18326 
.16628 
.15106 

.08746 
.08498 
.07883 
.06966 
.05833 

.15676 
.14173 
.13202 
.12782 
.12883 

.06334 
.04681 
.03097 
.01792 
.00937 

. 13818 
.12809 
.12110 
.11734 
.11676 

.04592 
.03358 
.02245 
.01357 
.00778 

.13422 
.14270 
.15269 
.16246 
.17037 

.00638 
.00919 
.01711 
.02861 
.04154 

.11911 
. 12397 
.13076 
.13879 
.14731 

.00557 
.00708 
.01205 
.01982 
.02941 

.17513 
.17587 
.17236 
.16501 
.15478 

.05351 
.06233 
.06639 
.06502 
.05861 

.15558 
.16286 

.16855 
.17216 
.17340 

.03967 
.04933 
.05722 
.06238 
.06417 

.14311 
.13164 
.12203 
.11563 
.11339 

.04856 
.03701 
.02643 
.01901 
.01635 

Integrals      -      Table   6 
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1 

=  2.5X 

t  = 

3.0\ 

t 

A 

B 

A 

B 

0 

. 33228 

-.31686 

.36421 

-.34776 

2 

.33339 

-.31569 

. 36569 

-.34622 

4 

.33669 

-.31212 

.37002 

-.34155 

6 

.34205 

-.30606 

.37705 

-.33356 

8 

.34929 

-.29731 

.38643 

-.32196 

10 

.35812 

-.28561 

.39773 

-.30638 

12 

.36813 

-.27072 

.41027 

-.28645 

14 

. 37882 

-.25237 

.42324 

-.26184 

16 

.38959 

-.23037 

.43560 

-.23234 

18 

.39967 

-.20462 

.44618 

-.19797 

20 

.40824 

-.17521 

.45359 

-.15914 

22 

.41436 

-.14245 

.45644 

-.11665 

24 

.41708 

-.10696 

.45341 

-.07185 

26 

.41550 

-.06968 

.44335 

-.02666 

28 

.40882 

-.03193 

.42555 

.01651 

50 

. 39649 

.00471 

.39992 

.05493 

32 

.37831 

.03835 

.36714 

.08578 

34 

.35453 

.06706 

.32883 

.10663 

36 

.32588 

.08900 

.28746 

.11585 

38 

.29369 

.10273 

.24624 

.11309 

40 

.25975 

.10744 

.20873 

.09957 

42 

.22627 

.10320 

.17836 

.07824 

44 

.19562 

.09111 

.15780 

.05345 

46 

.17003 

.07330 

.14834 

.03026 

48 

.15136 

.05275 

.14954 

.01351 

50 

.14067 

.03290 

.15906 

.00667 

52 

.13809 

.01713 

.17301 

.01085 

54 

.14264 

.00810 

.18670 

.02434 

56 

.15231 

.00724 

.19562 

.04294 

58 

.16438 

.01430 

.19658 

.06094 

60 

.17578 

.02739 

.18858 

.07284 

62 

.18369 

.04325 

.17318 

.07508 

64 

.18603 

.05797 

.15408 

.06735 

66 

.18195 

.06795 

.13609 

.05278 

68 

.17194 

.07080 

•12377 

.03695 

70 

.15781 

.06602 

.11994 

.02584 

72 

.14234 

.05520 

.12478 

.02362 

74 

.12856 

.04158 

.13569 

.03087 

76 

.11920 

.02923 

.14821 

.04434 

78 

.11595 

.02178 

.15744 

.05826 

80 

.11913 

.02139 

.15978 

.06664 

82 

A  L 

.12754 

.02806 

.15413 

.06598 

84 

.13883 

.03955 

.14235 

.05671 

86 

•a 

. 14996 

.05217 

.12855 

.04309 

88 

.15805 

.06184 

.11764 

.03128 

90 

. 16101 

.06545 

.11351 

.02665 

50 

APPENDIX  C 

The  Fresnel  Integral  Solution 

This  approach  makes  use  of  the  integral  equation  described  at  the 
beginning  of  Appendix  A 


2  -  k  2  )  f 

—~—  J  E(Q)  E(Q2\kJ))    da 


(k  "    -   k 
,1. 


For  the  iteration  series  solution,  the  field  at  a  distant  point  P  was 
calculated,  using  the  infinite  line  source  at  the  origin  as  the  primarjr 
source.   For  the  following  calculations  the  principle  of  reciprocity  will 
be  employed  merely  to  illustrate  that  there  exist  two  equivalent  view- 
points from  which  the  same  radiation  pattern  formulation  may  be  obtained. 
For  the  reciprocal  case,  the  quantity  to  be  obtained  becomes  the  elec- 
tric field  at  the  origin  due  to  an  incident  plane  wave.   In  other  words 
the  point  P  will  be  taken  to  be  the  origin  and  the  field  in  the  dielec- 
tric sheet  [E(Q)]  will  be  due  to  the  incident  plane  wave.   The  radiation 
pattern  is  thus  given  by  the  absolute  value  of  E(P) . 

As  in  Appendix  A,  it  will  be  assumed  that  the  dielectric  sheet  is 
very  thin  so  that  the  integral  over  its  cross-section  can  be  replaced 
by  its  thickness  multiplied  by  the  line  integral  along  the  sheet.   Further- 
more it  will  be  assumed  (also  as  in  Appendix  A)  that  the  field  inside 
the  sheet  is  the  same  as  the  incident  field  in  free  space.   The  incident 
field  E"  is  a  plane  wave  with  arbitrary  phase  and  magnitude.   It  is  con- 
venient to  choose  the  phase  and  magnitude  so  that  the  field  at  the  origin 


51 


has  unit  magnitude  and  zero  phase  angle,  that  is  E  (P)  =  1.   This  gives 

jk  x'cos  0 
o 
the  electric  field  along  the  sheet:   E(Q)  =  e  .It  must  also 

be  observed  that  p   (the  distance  between  P  and  Q)  is  now  |x|  so  .that 


,(2) 


,(2), 


H  2  (k  p)  becomes  IT2^(k  Ix'l). 
o    o  o    o  ■ 

The  notation  in  the  remainder  of  the  derivation  will  be  made  slightly 
more  compact  if  all  distances  are  in  free-space  wavelengths  (indicated 

—   x  —  i 

by  a   bar   over   the   distance   symbol).      That    is   x  =  r,    i    =  r, 

eJk0x'cos  <p=   eJ2li'coS  <p     ^  H(2)  ,       =  H(2)(2ff|-lK 

o     o  I   '      o 

Now  we  can  write  the  integral  relation  giving  the  field  at  the  origin 
P  due  to  an  incident  plane  wave. 


E(P)  =  1  -  JTTCe   -  1)  t 
r 


i 

ej27rx^cos^H(2)(27r|-/|) 


dx' 


-i 


The  primes  (previously  used  to  indicate  coordinates  in  the  dielectric 
sheet)  can  now  be  dropped  since  there  are  no  unprimed  coordinates  in  the 
formula. 

So  far  the  procedure  has  been  essentially  equivalent  to  that  used 
in  the  iteration  series  approach.   Now,  in  order  to  make  use  of  the  Fresnel 
integrals,  it  is  necessary  to  write  the  above  equation  in  the  form 


E(P) 


jff2(€ 


-  1)  t 


eJ27Tx  cos  %(2)(2ffj-()  d- 
o 


-CO 


r 


o  J 


-OO 
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This  can  be  written  as  E(P)  =  1  +  E  (P)  -  E  (P)  where  E  (P)  is  the  field 

12  1 

due  to  an  infinite  dielectric  sheet  and  E  (P)  is  due  to  a  pair  of  semi- 

2 

infinite  sheets  (one  from  -i  to  -°o  and  the  other  from  +i  to  -h»)  . 

Calculation  of  E  (P) 

_ 1 

If  we   set   (3  =   -   jlT2(6      -   1)    t,    we   have 


Ei(P)    =   (3     J       H^)(27T|x|)    ej27rX   C°S   ^  dx 

-OO 


=    2(3  H^2    (27Tx)    cos(27Tx  cos   (p)    dx 


Setting  u   =   27Tx,    we  can  write 

00 

E   (P)    =  |  H(2)(u)    cos(u   cos  (p)    du 

1  *  J         ° 

o 

It  is  worth  noting  that  the  range  of  integration  is  from  zero  to  infinity, 
In  the  iteration  series  derivation^  the  integration  did  not  start  from 
zero  but  from  a  small  value  6.   This  took  into  account  the  fact  that  the 
experimental  dipole  was  inserted  in  a  hole  of  radius  6.   For  the  deriva- 
tion to  follow,  the  hole  will  not  be  considered  (6=0)  in  order  to  take 
advantage  of  the  following  integral  formulae: 
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oo 
|        J    (ax)    cos(bx) 


dx  =  — 


Jo  fa2    -  b2 


0  <  b  <  a 


=   0 


b  =  a 
a  <  b  <  oo 


OO 


J 


N   (ax)    cos(bx)    dx  =   0 
o 


0  <  b  <  a 


=  oo 


b  =  a 


v,2  o2 

b     -   a 


Using    the  above   formulae,    we  can  write 


a  <  b  <  oo 


V«  =1 


/l   -   cos2    (p 


W   sin  <p 


-    j7T(€      -    1)    t 
r 

"sin  (p 


which  is   the   required   expression   for   the   scattered   field   due   to  an   infinite 
dielectric   sheet. 


Calculation  of  E 


Setting  S  =   -   jIT   (£      -   1)    t,    we  have 

r 


.oo 


E   (P) 

2 


H(2)(2ir|x|)    eJ2?rX  C°S  *  dx  +  f     H(2)(21Tix 
o  >  o 


j21Tx  cos  <p    — 
e  dx 


-oo 


f 


.OO 


I    I  i^2)/07r~\       -j2fx  cos   (p    -       J  1I(2)l-«--x      J2JTx  cos   <p    - 

=   p  <  H        (23Tx)    e                             dx  +  H        (27Tx)    e                           dx 

■>_        °                                                             >_  ° 

i  i 
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If  1    is  about  one-half  or  greater  (as  in  the  experiments  performed),  the 
Hankel  function  can  be  very  closely  approximated  by  the  following  formula: 


H(2>(27Tx) 
o 


-  E 

1     e-j(27Tx  -  4) 


7T  1/  x 


Making  this  substitution,  we  have 


!T    f     ,oo 


E    (P) 


6      J    4 

S   e 
IT 


-j27Tx(l+cos  (p) 


OO 


-j27Tx(l-cos   (p) 


dx  + 


/I 


j  ?rr°°  ~JV 


,oo     -jy  x 


6          4)   |       e                          |       e 
-  «        1  dx  +         dx 


7T 


i  /x 


/3 


where  -y     =   27T(1    +  cos  <p)    and  ^     =   27T(1    -   cos  <£>)  . 

It   is   now  necessary   to   evaluate   integrals   of    the    form 


.00  .  JX5 

£ dx  = 


/x 


cos   vx 


^ 


dx   -   j 


si 


IL^S 


dx 


^  i/x 

a  v 


These  will    be  expressed   in   terms    of  the   Fresnel    integrals  C   and  S   which 
are  defined  as    follows    (see   reference   7): 


•<f 


1 


cos   u     ,        -.  f   /2v 

du,    S(r¥ 

u 


sin  u 


JZR J  ./~lT 


du , 
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Since   S(°oj    =  C(°o)   =   2,    it    is   possible   to  change   the   range   of    integration 
to  obtain 


cos   u 


/T 


du   =      \[W 


2        C  W7T   )\ 


/} 


I  fr 


sin   u     , 
— —  du   = 


2  ~  s(/lf 


Returning   to   the   integral    form   found   in   the   problem,    we  make   the   transfor- 
mation -yx  =   u  which  gives 


e~^X  1 

- —   dx  = 


-ju 


/x"  A 


du  = 


\a 


r 


\a 


00 


00 


cos    u     ,          .       1 
— du   -j   


sin  u 


du 


Va 


277 


J  ; 


277 


V7 


I  _  c  ( ,^a 

2         U  U     77 


F   (\a)    -j    F   (\a) 
c  s 


where 


Fc(\a) 


r  f    2Ya 


F   (7a) 
s 


-  s 


\  V    77 


The   functions  C    and   S   are    tabulated   in   reference   7    for  y&   from  0   to   50 
in    ,01    steps , 

It    is   now   possible   to   obtain   the   field  E    (P)    in   terms   of   the   tabi- 

2 

lated   functions. 
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E    (P) 
2 


1     f 
7T 


F    (y    i) 
C     «l 


J    F   U  i) 
s     »1 


f27T 


F    (v  4) 
C      »2 


j    F    (v  i) 
S     "2 


j"  I  F    (v  i)      -    j    F    (v   i)  F    (v   |)    -   j    P    (v   i) 

3     „  4  J     c    Ti s    'i  c    *2 s    '2 

_  ^      e      >  ______ — , _ . +  _________________ 

/l    +  cos  0  \/l    -   cos  (p 


$ 


■    =^-T ^     1    l/l     ~     COS     #> 

JT   sin  <P 


F   (v  i)    -j    F   (v  i) 

c    'i  s    »_ 


+        /l     +    COS    <£> 


F   (v  i)    -jF   (v  i) 
c     l2  S     J2 


Now      ,/l    +   cos  #>  =      /2 
J 

J4        1    +   j 
^       _  — . — _y.^      Thus   we  have 


<P 


<P 


cos  ^  and      j/1    -   cos   (p  -      /_~  sin  ~  and   in  additi 
_*  — 


on 


/¥ 


,  (p,  .  g_M_>    (sin  | 

2  TT   sin  <£>  j  2 


F   (v  i)    -jF   (vi)      4    cos  * 
c    li  s    •_  2 


^ 


F    (v   i)    -jF    (v   i) 
C     '2  S     '2 


=  |^^(sinf  F(Yi)    +cosf  F(VJ) 
7T   sin  </M  2      c    *i  2      c    T2 


<P 


<P 


sin  £  F   (v  i)    +  cos  ~   F   (v  i) 
2      s    'i  2      s    '2 


=— ?— -j,]  sin  f  Ff   (v  4)    ■  (v  i) 

AT   sin  <jP  2  [  c    Ti 


s    «i 


+  cos 


0 


F    (v  i)    +   F   (v  i) 
c     \2  S     '2 


+  j    sin  - 


<a 


F   (v  i)    -    F   (v  i)      +  j    cos 


<P 


C      »1  S      Jl 


F    (v  I)    -    F    (v  i) 

C     ?2  s     "2 


=   S — ? — 7n  i  N   +   jM  [      =    " 
IT    sin  &  I 


-^h1^    JN    +  JM 

sin  <p  I 
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where 


N   =    sin 


<Pr 


F    (Y.|)    +    FCy.l) 
c    '1  s    '1 


sin  - 


1    -sW-y 


/-  c  v 


2Yl! 

77~ 


+    cos  - 


+    cos  I 


Fc(v2'i)  +  fs(y2o 


/ 


1    -    S 


2Vr 


77 


'-  cv 


^7i 


77       / 


M  =    sin  — 

2 


Fc(Yli)    "    Fs(Yli) 


+     C°Sf    IVV*     "     FB«V>] 


sin  ~ 
2 


(S)- 


2Yi* 


+    cos   ~- 


PV     -  [  /V  \ 
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77 


y   J 


<y      =    27T(1    +    COS    Cp) 
Y2   -   277(1    -    cos   <p) 

Thus,  with  the  help  of  a  table  of  Fresnel  integrals  S  and  C,  the  values 
of  M  and  N  can  be  calculated  for  any  given  0  and  (p„ 

The  Total  Field  E(P) 

E(P)  is  given  by  the  relation 


=  1  +  —  •    .      -  £1  •  — : 1  N  + 


E(P)  =  1  +  E.,  (P)  -  E  (P)  =  1  +  £  •    .     -  £  •  — : -,  N  +  JM 

1       2  7T    sin  (p       77    sin  (£> 


s  ir(£  -  i)t 

:  X  +  FTT^  tl  ■  N  -  J«]  n  1  -    sr±n(p  [M  +  j(l  -  N>] 


The  radiation  pattern  P(cp)    is  given  by  the  absolute  value  of  E(P') 
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APPENDIX  D 

The  Asymptotic  Form  of  the  Fresnel  Integral  Solution 

One  form  of  the  Fresnel  integral  can  be  expressed  asymptotically  as 
follows: 


oo  o 


jb    /27T 


1    - 


2b' 


-Jb2 


jb  fzn 


if   2b*   »  1 


2  2 

Making  the  transformation  t  =  "Yx  and  b  =  Va,  we  obtain 


-JVx       -jYa 


dx 


a  {x  J7  /a" 

if  2va  »  1 

From  Appendix  C  we  have  an  expression  for  E  (P),  the  scattered  field  at 

2 

the  origin  when  a  plane  wave  is  incident  on  two  co-planar  semi-infinite 
dielectric  sheets. 


7T 


_£   '4 


'  J»     -j"Y  x 


E  (P)  =Je 
2        « 


dx  + 


oo  -jY  x 
e   2 


.  i     /x 


dx 


/I 


where  v     =   27T(1    +  cos   (p) ,    v     =   21T(1    -   cos  (p)  .      If   2v  i    »  1    and  2v  i   »  1 
ll  '     T2  1  2 

we  can  write 
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E   (P)    =  |   e   4  < 
2  " 


-JY*         -JY  J 

e       1  e       2 

4    


JY,  M      jyo  M 


1T(€      -   l)t        J4       f    -j27Ti(l+cos  (?)  -j21Ti(l-cos   <p) 


-nil 


1    +   cos    </> 


1   -   cos  cp 


—        U" 
H(€r    -    l)t      e"j(2iri    '   4)     [     -jOTi    cos   <p        e+j27Ti    cos   <? 

1      1    +  cos   <p       +      1    -   cos  (p 


irvJl 


(o)  - 

Now   the   large-argument   approximation   to   the  Hankel    function  H        (21Ti)    is 
given  by 


H(2)(23Ti) 
o 


■j(27Ti    -  H) 


=   R(i) 


ir  i/| 


Thus   we   can  write 


7T(€      -   i)t 

E    (P)    = £- R(i) 

2  2 


-j27Ti    cos   <£>  +J27T!    cos   <p 

1    +  cos   (p  1    -   cos   <£> 


It  has  already  been  shown  in  Appendix  A  that  E  (P)  is  given  by 

1 


7T(€   -  l)t 

E  (P)  =  -  j  7 -r- - 

I  sin  (p 


The  total  field  at  the  origin  E(P)  now  can  be  obtained  using  the  asymp- 
totic formulae.   We  have 
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E«,P)  =  1  +  E  (P)  -  E  (P) 
1        2 


=  1  + 


He     .   1;  t"  Lj_  +  SO)  \e^   -  *  +  ggg  -  t\ 

r         sin  (p         2   [  1  +  cos  <P      1  "  c°s  <P  J 


This  formula  has  been  derived  using  the  asymptotic  expansion  for  the 

Fresnel  integral  and  holds  if  2y  i   »   1  and  2-y  i  »  1.   That  is 

1  2 

47Ti(l  +  cos  (p)   »   1  which  means  that,  for  a  fixed  value  of  S. ,    <P   is  not 

o       o 
permitted  to  approach  0   or  180  . 

-1 
In  the  above  formula  the  factor  (j  sin  (p)        comes  from  an  integration 

over  an  infinite  dielectric  sheet.   The  other  term 


R(j) 

2 


~j2Mi    cos  (p         +j2ffi    cos   <p 

_   + 

1   +  cos  <p  1    -   cos  cp 


suggests  that  the  thin  dielectric  strip  may  be  represented  by  two  infinite 

line  sources  in  free  space.   These  equivalent  line  sources  are  located 

at  x  =  -i    and  x  =  +i .   Furthermore  these  equivalent  line  sources  have 

individual  radiation  patterns  of  the  form  (1  +  cos  (p) 

Now  E(P)  is  of  the  form  E(P)  =  1  +  (€   -  1)  t  G  where  G  is  a  complex 

r 

number  and  a  function  of  <P  and  I.   If  (€   -  1)  t  is  very  small  so  that  the 

r 

magnitude  of  (€   -  1)  t  G  is  much  smaller  than  one,  then  the  radiation 
r 

pattern  given  by  P(<#)  =  il  f  (€   -  1)  t  Gl  can  be  approximated  by 

P(<p)  =  1  +  R  \  (€   -  1)  T  0  j  where  R   indicates  "the  real  part  of".   Thus 

the  radiation  pattern  formula  becomes 

P(<0  =  1  +  ?T(€   -  1)  t  R  )  ~^-  |  ~ -r-    +  ^- rg-    } 

r         e  (  2   L.  1  +  cos  cp  1  -  cos  <p    J 
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The  expression  inside  the  curly  brackets  is  the  line  source  term  described 
in  the  previous  paragraph.   Thus  a  thin  dielectric  strip  with  a  low  dielec- 
tric constant  can  be  represented  approximately  by  two  line  sources  in  free 
space . 

The  "Two-Source  Equivalent"  with  Directivity 

This  report  deals  with  values  of  i  which  are  integral  multiples  of 
a  half  wavelength.   Thus  we  may  set  St    =  —  where  n  is  an  integer.   The  above 
formula  for  V(Cp)    can  be  simplified  and  written  as 


.-  r |  cos(n7T  cos  Cp)    -  cos  cp   sin(n7T  cos  Cp) 


/n"  sin2  Cp 


where  the  directivity  factors  (1  +  cos  cp)        have  been  included. 
The  "Two-Source  Equivalent"  without  Directivity 

If  cp   is  further  restricted  to  be  not  far  from  90  }    then  we  may  make 

the  approximations  sin  <p  =  1  and  cos  Cp  =  0 .   This  means  that  the  direc- 

-1 

tivity  factors  (1  +  cos  Cp)        become  unity  and  the  dielectric  strip  is 

represented  by  two  ordinary  non-directive  line  sources.   The  radiation 
pattern  formula  becomes 


-  (-Dn 
V{cp)    =  1  +  (€   -  1)  t —  cos(n7T  cos  Cp) 

•/n 


where,  as  before,  n  is  an  integer  and  gives  the  number  of  half -wavelengths 
in  the  length  St  . 
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APPENDIX  E 

When  an  electromagnetic  wave  propagates  along  a  dielectric  slab,  the 

wavelength  is  shorter  than  in  free  space.   The  ratio  of  guided  wavelength 

X  to  free  space  wavelength  X  can  be  obtained  by  the  source-free  solutions 
g 

of  Maxwell's  equations  for  an  infinite  dielectric  slab  in  space.   Such 
solutions  result  in  the  following  characteristic  equations  for  the  TE  modes: 


r  kxti 
tan  L  2  J 

/k  2    -   k  2 

/     1              ° 

1     ** 

1 

cot 

'k   t" 

X 

2  _ 

/k  2    -   k  2 

/    1             o 

"V      K* 

-    1 

The  propagation  constant  along  the  dielectric  slab  is  given  by 


I    V  1      x 

and  the  wavelength  ratio  is 


\    k 
_J>  =  _2 
\  k„ 


where 

t  =  dielectric  slab  thickness 
27T 


k  =  w 
o 


/m  €   =  -j-  -   the  propagation  constant  in  free  space 

V  "O  O     A. 


{*)  /jx  €  €   a  the  propagation  constant  in  a  dielectric 
medium  of  relative  dielectric  constant  € 


r 
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27T 
k  =  -j—   =  longitudinal  propagation  constant  for  a  dielec- 

g 

trie  slab  in  space. 


For  very  thin  dielectric  slabs,  only  the  lowest  order  TE  mode  can 

propagate  and  this  mode  has  no  cutoff  frequency.   For  the  lowest  order 

\ 

or 

TE  mode,  the  propagation  constant  k   (and  consequently  the  ratio  -r^)  can 

x  At 

be  obtained  from  the  solution  of  the  transcendental  equation 


tan 


k  t 
x 


k  2  -  k  2 


k  2 

X 


\ 


cr 

For  the  case  of  a  polystyrene  slab  (€   =  2.55),  the  values  of  -r—   are  tabu- 
lated below. 


t 

t 

(in  inches 

g 

\ 

at 

5000  mc) 

1/4  " 

\ 

.1058 

.908 

.05292 

1/8  " 

.972 

.02646 

1/16" 

.994 

.01323 

1/32" 

.998 

If  t  is  assumed  to  be  very  small,  the  above  transcendental  equation 
can  be  put  into  a  very  simple  form.   The  use  of  a  trigonometric  fundamental 
identity  changes  the  equation  to 


k  t      k  2 
2   x        x 
cos 


2    k  2  -  k  2 


The  half-angle  formula  gives 
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2k 


cos  k  t  = 


x     ,2    ,2 
k   -  k 
1      o 


-  1 


The  assumption  that  t  is  small  permits  the  approximation 


2    9. 

k  *    V" 


cos  k  t  =  1  - 
x 


which,  puts  the  equation  into  the  form 


\  k  ' 
\  1 


i         t2  v1 


Now  k 


k   -  k„   and  k 


2       2 
k  '  =  k   (€   -  1) .   Making  these  substitu- 
o     o   r 


X  l  JL  1 

tions  and  appropriate  binomial  approximations  (assuming  t  to  be  small) 


we  get.  the  approximate  expression  for  k  , 


k„  =  k 
i     o 


!  ,£!((£   .  1}  ".2 


2  L   r 


* 


from  which  the  desired  values  of  -r2  can  be  calculated.  These  are  given 
in  the  following  table  and  compared  with  the  results  of  the  direct  com- 
putation of  the  transcendental  equation  solution,   (€   =  2,55). 


t 

Exact 
.908 

App 

roximate 

.1058 

.868 

.05292 

.972 

.967 

.02646 

.994 

.992 

.01323 

.998 

.998 
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There  are  two  interesting  similarities  between  the  above  approximate 

formula  and  the  radiation  pattern  formula  (iteration  series  solution). 

Both  formulae  are  accurate  over  the  same  range  of  €   and  t,  and  the  factor 

r      ' 

(€   -  1)  y  is  prominent  in  both  expressions, 
r      a. 


